Exact solutions via Prabhakar fractional approach to investigate heat transfer and flow features of hybrid nanofluid subject to shape and slip effects

The core devotion of this study is to develop a generalized model by means of a recently proposed fractional technique in order to anticipate the enhancement in the thermal efficiency of engine oil because of the dispersion of graphene and magnesia nanoparticles. In addition to investigating the synergistic attributes of the foregoing particles, this work evaluates shape impacts for column, brick, tetrahedron, blade, and lamina-like shapes. In the primary model, the flow equation is coupled with concentration and energy functions. This classical system is transmuted into a fractional environment by generalizing mathematical expressions of thermal and diffusion fluxes by virtue of the Prabhakar fractional operator. In this study, ramped flow and temperature slip conditions are simultaneously applied for the first time to examine the behavior of a hybrid nanofluid. The mathematical analysis of this problem involves the incorporation of dimension-independent parameters into the model and the execution of the Laplace transform for the consequent equations. By doing so, exact solutions are derived in the form of Mittag–Leffler functions. Multiple illustrations are developed by dint of exact solutions to chew over all aspects of temperature variations and flow dynamics. For the preparation of these illustrations, the details of parametric ranges are as follows: \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.00 \le \varUpsilon \le 0.04$$\end{document}0.00≤Υ≤0.04, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2.0 \le Gr_1 \le 8.0$$\end{document}2.0≤Gr1≤8.0, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.5 \le Sc \le 2.0$$\end{document}0.5≤Sc≤2.0, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.1 \le \uptau \le 4.0$$\end{document}0.1≤τ≤4.0, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.1 \le d \le 0.6$$\end{document}0.1≤d≤0.6, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.2 \le \lambda _1 \le 1.5$$\end{document}0.2≤λ1≤1.5, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.0 \le Gr_2 \le 7.0$$\end{document}1.0≤Gr2≤7.0. The contribution of differently shaped nanoparticles, volume proportions, and fractional parameters in boosting the heat-transferring attributes of engine oil is also anticipated. In this regard, results for Nusselt number are provided in tabular form. Additionally, a brief analysis of shear stress is carried out for fractional parameters and various combinations of magnesia, graphene, and engine oil. This investigation anticipates that engine oil’s hybridization with magnesia and graphene would result in a 33% increase in its thermal performance, which evidently improves its industrial significance. The enhancement in Schmidt number yields an improvement in the mass transfer rate. An increment in collective volume fraction leads to raising the profile of the thermal field. However, the velocity indicates a decreasing behavior. Nusselt number reaches its highest value (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Nu=8.1363$$\end{document}Nu=8.1363) for the lamina shape of considered particles. When the intensity of the buoyancy force augments, it causes the velocity to increase.

www.nature.com/scientificreports/ Magnesia and graphene nanoparticles have certain significant features that make them useful for a variety of engineering and industrial applications. For instance, graphene possesses a large surface area to volume ratio, which makes it ideal for use in applications where surface area is important, such as catalysis and heat transfer. Graphene nanoparticles are highly conductive, meaning they can easily conduct heat and electricity. Due to this feature, they are highly effective for electronics and energy storage applications. These particles are incredibly strong and rigid therefore, their use for reinforcement in coatings and composites provides desirable outcomes. Furthermore, they can be utilized in sensors, drug delivery, water filtration, and tissue engineering. On the other hand, magnesia is thermally stable at extremely high temperatures because it offers significant resistance to current and has a high potential for heat conduction. These particles are also biocompatible therefore, they are effective for certain biomedical applications like imaging and cancer therapy. Also, the production of specific optical materials that are used to treat indigestion and heartburn involves magnesia. Since magnesia offers resistance to moisture and fire, it is one of the fundamental components of construction materials.
The concept of adopting fractional methods to modify regular models gave rise to a new field called fractional calculus. Due to its extensive implementations in diverse real-life conditions, it is a discipline that is now undergoing rapid evolution. In recent times, multiple experts from various scientific fields have reported that results procured by the virtue of fractional approaches are more reliable, and the execution of fractional operators for modeling purposes ensures the specificity and accuracy of outcomes. In addition, they offer a more precise interpretation of the under-observation process. Besides, a fitting adjustment of fractional parameters leads to a good accordance between theoretically derived solutions and experimentally established findings. These additional benefits have motivated a number of researchers to examine physical problems in fractional environments and perform comparative investigations. The usages of fractional models can be found in a variety of areas such as dynamical systems, control theory, disease and population modeling, electromagnetics, economics, fluid mechanics, mathematical biology, and so forth. Regarding fluid mechanics, memory and self-similar effects of fractional derivatives are significantly crucial to fully comprehend the rheological features, thermal performances, and viscoelastic behaviors of fluids. So far, a variety of fractional operators composed of various mathematical formulations have been presented. Each of them has unique limitations and benefits. In this list, Caputo and Riemann-Liouville are the most frequently utilized operators, and their formulations involve a power-law kernel 32 . The other well-known operators are Prabhakar, Atangana-Baleanu, Hilfer, Caputo-Fabrizio, Hadamard, and Grüünwald-Letnikov fractional operators [33][34][35] . As compared to standard methods, modern-day researchers are preferring fractional modeling techniques in order to provide more authentic descriptions of physical mechanisms based on generalized solutions. Several fractional operators have been utilized until now to chew over the complexities of diverse natural phenomena. Fallahgoul et al. 36 analyzed the impacts of vital characteristics of fractional processes, for instance, self-similarity, path-dependency, and long-range memory, on financial theory, economics, and financial models. Sinan et al. 37 used the Atangana-Baleanu operator to establish a model for an in-depth investigation of the malaria disease. They discussed the effectiveness of precautionary measures and medication for the reduction of the disease's spread. Asjad et al. 38 evaluated the control of generalized boundary conditions on heat-conducting properties of nanofluids by means of a fractional system. Raza et al. 39 explained the thermal aspects of water and kerosene oil based nanofluids with the aid of semi-analytic fractional solutions. Ikram et al. 40 established multiple fractional models to explore the heat-transmitting efficiency of hybrid nanofluids during channel flows. Some of the latest studies in which thermal and flow distributions are examined via fractional operators can be viewed in [41][42][43][44] .
A meticulous analysis of the literature reveals that there are not enough studies on such hybrid nanofluids, which contain oils as host fluids. This research gap further expands if the utilization of fractional operators and the derivation of exact solutions are simultaneously taken into consideration. In addition, it is noted that shape aspects of nanoparticles have not been given adequate importance because most of the reported investigations communicate results for spherical structures. This work is an attempt to address all these concerns. The principal feature of this analysis is to investigate the consequences of hybridizing magnesia and graphene nanoparticles with engine oil. The shapes' influences are given significant attention. In this regard, it is assumed that the observed particles have blade, column, lamina, brick, and tetrahedron-like shapes. This research work also aims to explain the flow dynamics and thermal behavior of the consequent hybrid nanofluid in terms of a fractional model. To serve this purpose, generalized relations for thermal and diffusion fluxes are established with the aid of the Prabhakar fractional operator. The inclusion of dimension-independent quantities in the principal system lays the foundation for the implementation of the fractional operator. In this work, the behavior of hybrid nanofluid for simultaneous application of ramped flow and thermal slip conditions is examined for the first time. To solve the ensuing fractional system, the Laplace transformation is executed, and exact solutions are produced in the form of multi-parametric Mittag-Leffler functions. Various tables and graphical illustrations are presented to effectively evaluate flow patterns, thermal profiles, shape impacts, contribution of influential parameters, concentration field, and heat transfer performance. The illustrations are compared for slip and no-slip cases and for lower and higher time values to emphasize the importance of slip conditions, transient effects, and the ramped velocity condition. Some modifications in fractional parameters are performed to retrieve thermal and flow functions for the classical case, and their graphical comparison is carried out with those acquired via the fractional model.

Statement and model formulation
In this work, the hybridization of multi-shaped magnesia and graphene nanoparticles with engine oil is examined to discuss the significant modifications in thermal and flow characteristics of engine oil. The geometrical arrangement for this analysis involves an infinite upright wall that serves as the solid-fluid interface. Initially, the temperature of hybrid nanofluid is ∞ with a uniform concentration C ∞ , and the system demonstrates no motion. Later on, ramped movement of the bounding wall and a temperature variation due to slip effects www.nature.com/scientificreports/ disturb the system. In the mathematical sense, a piece-wise function describes the ramped movement in such a way that the velocity is a time-dependent function ( U 0 ( τ/τ 0 ) ) for a particular duration ( τ ≤ τ 0 ) . After that (for τ > τ 0 ), the velocity has a constant value ( U 0 ). Meanwhile, the concentration changes from C ∞ to C w . Far from the wall, the flow function associated with hybrid nanofluid attains a zero value, and thermal and concentration functions again achieve ambient values ( ∞ and C ∞ ). Figure 1 provides the geometrical setting of this study. The mathematical model is developed considering the following assumptions • The flow, concentration, and thermal functions only contain one axial component ( Ψ ) due to the infinite length of the bounding wall. • The flow is one-dimensional and unidirectional.
• Engine oil is in thermal equilibrium with magnesia and graphene nanoparticles.
• Nanoparticles are assumed to have brick, column, lamina, blade, and tetrahedron shapes.
• The viscosity dissipation does not disturb the heat transfer process.
• The buoyancy effects are addressed using the Boussinesq approximation 45 .
Taking the aforestated assumptions and description into consideration, the principal equations for this problem are derived as 46 Fourier's law for thermal flux and Fick's law for diffusion equation are respectively provided as  www.nature.com/scientificreports/ The following conditions are associated to the above-mentioned system of equations

Mathematical relations for thermal and physical properties
The major factors that make a hybrid nanofluid preferable to a standard industrial fluid are its enhanced physical and thermal characteristics. These features of the involved nanoparticles significantly influence the development of flow and have a substantial impact regarding the thermal usability of the consequent hybrid nanofluid. However, when it comes to describing all these features mathematically, there is no universal model that can address them simultaneously. However, a number of researchers tried to study different aspects of these features through experiments. As a result, several mathematical relations have been established to explain the thermophysical attributes of nanoparticles. Later on, to adequately characterize these properties of hybrid nanofluids, these models were effectively adapted. This section primarily focuses on outlining a fundamental mathematical relationship for each thermo-physical feature and its alteration in the case of hybrid nanoparticles.
Viscosity. Various factors, such as density, drag effect, and viscous forces, contribute to the formation of flow patterns. Additionally, it is crucial to anticipate the ability of a nanofluid to oppose the deformation. In this regard, Brinkman proposed a model in 1952 47 that received significant attention later, and it is the most frequently applied model in recent times. It anticipates the viscosity as The expression for the viscosity of the considered hybrid nanofluid is defined as Specific heat capacity and thermal and concentration expansion coefficients. The following equations are used to explain specific heat capacity and thermal and concentration expansion coefficients The altered forms of the above-mentioned relations for hybrid nanofluid are as follows 10 1 − Υ Gra 25 10 . www.nature.com/scientificreports/ Density. The expression for an adequate estimation of the density is provided as The modification of the above-mentioned expression for magnesia and graphene based hybrid nanofluid is presented as Thermal conductivity. In general, thermal conductivity has a strong impact on the efficiency of a nanofluid. Thus, accurate measurement of thermal conductivity is quite essential. The model introduced by Hamilton and Crosser 48 is currently the leading model for the quantification of thermal conductivity. This model also counters shape influences effectively therefore, most of the researchers prefer it when objectives include the evaluation of shape effects. This model relates the thermal conductivity of nanofluid, host fluid, and nanoparticles, and the shape factor in the following way where the choice of a specific shape of nanoparticles determines the value of the shape factor (s). The working fluid in this study is composed of hybrid particles therefore, the extended version of Eq. (19) contains two parameters s 1 and s 2 , which correspond to graphene and magnesia nanoparticles, respectively. The extended version is communicated as where In this section, the subscript "hnf" stands for hybrid nanofluid, "nf" denotes nanofluid, and "hf" symbolizes host fluid. For computation purposes, the values of s 1 and s 2 will be chosen from Table 1. The computational values for thermo-physical features are accessible from Table 2.

Exact solutions of a non-dimensional generalized model
In this section, firstly, the developed model will be made dimension-independent to provide the basis for the application of the fractional derivative. This purpose will be achieved by plugging some unit-independent parameters into fundamental equations and connected constraints. Secondly, expressions for thermal and diffusion fluxes composed of generalized Fourier's and Fick's laws will be incorporated in the subsequent dimensionindependent system to procure a fractional model. Lastly, a comprehensive mathematical analysis will be carried (15)  www.nature.com/scientificreports/ out for the computation of exact solutions. In this analysis, the fractional model and connected constraints will be treated with the Laplace transform (LT). To achieve the first goal, new quantities are presented as follows The substitution of thermo-physical expressions and the above-mentioned quantities in Eqs. (1)-(5) returns The unit-free forms of pertinent conditions are imparted as where Table 3 provides the parameters that appear in Eqs. (23)- (27). Now, this dimension-independent classical model will be shifted into a fractional setting by generalizing the expressions of thermal and mass fluxes using the Prabhakar fractional derivative. These expressions are receptively supplied as where the Prabhakar fractional operator ( D η ζ ,σ ,ε ) for an arbitrary function G(u) is given as 34 www.nature.com/scientificreports/ where is the Prabhakar integral. The Mittag-Leffler function composed of three parameters and the Prabhakar kernel are given as 51 The application of LT on the Prabhakar fractional operator yields the following expression Mathematical analysis of temperature function. The new forms of temperature equation (Eq. (24)), Fourier law (Eq. (31)), and connected constraints (Eqs. (29) and (30)) that derived through the utilization of LT are respectively imparted as Quantities relating thermophysical features and volume fraction of nanoparticles The Laplace domain concentration function is evaluated considering the relevant conditions (Eq. (45)), and it is supplied as  (48) (52) www.nature.com/scientificreports/ , www.nature.com/scientificreports/

Results and discussion
The principal purpose is the assessment of enhancements in the thermal properties of engine oil as a result of the immersion of magnesia and graphene nanoparticles. The Prabhakar fractional operator is exercised as a generalization tool to establish fractional versions of classical equations. Uniform concentration and slip temperature conditions are considered together with the ramped motion of an infinite vertical bounding wall. The ramped velocity function and buoyancy forces (mass and thermal) are the major factors leading to the instigation of the flow. The governing system of equations is composed of flow, concentration, and energy functions, and it is a partially coupled system. Utilizing the Laplace transform, the generalized system is solved, and solutions comprised of Mittag-Leffler functions are derived. This section is organized to provide these solutions in tabular and graphical forms, which were obtained using MATLAB. The graphical illustrations are presented for slip and no-slip cases and for lower and higher time values. A comparison of classical and fractional solutions is also communicated graphically. Additionally, Nusselt number and skin friction coefficient are comprehensively investigated to analyze several phenomena, for instance, the impacts of column, brick, tetrahedron, blade, and lamina shaped particles on thermal efficacy, escalation in heat transfer rate, and the consequence of altering fractional parameters and volume proportions on shear stress. The focus of Fig. 2 is to discuss the implications of altering the fractional parameter σ . Figure 2a,b indicate a substantial decline in the outcomes of thermal and concentration functions in response to minor enhancements in σ . For the thermal field, the consequences of modifying σ are same whether slip effects are taken into consideration or not. However, when the slip condition is applied, the corresponding temperature graph is always lower than the one associated with the temperature function obtained without the slip condition. The aforesaid finding is likewise valid for velocity graphs. Moreover, an interesting behavior of the flow function is noticed on analyzing the involvement of the parameter σ . The variations in the flow field for time-dependent and uniform conditions are opposite to each other. The flow profile rises for the ramped case, but subject to the constant condition, the velocity field displays dropping profiles, as demonstrated in Fig. 2c,d. The substantial disturbances in flow, concentration, and heat functions against slight changes in σ show that fractional operators are highly efficient for controlling these functions according to the physical situations. Furthermore, the www.nature.com/scientificreports/ parameter-adjusting property of these operators ensures the acquisition of agreement between theoretical and empirical results. So far, an analysis is made on the basis of a single parameter. However, investigating the behavior of principal functions for modifications in all the involved fractional parameters at the same time is equally essential. For the current study, this task is achieved by plotting the graphs of flow, concentration, and heat functions in Fig. 3 for various values of parameters ζ , σ , and η . From Fig. 3a, it is discovered that an increase in the aforementioned parameters reduces the heat function. The response of the concentration field is also the same, as depicted in Fig. 3b. Figure 3c,d disclose that the velocity field behaves in an opposite manner for uniform and time-dependent conditions. It reduces for the earlier case and increases for the later case. The three-parametric kernel, which facilitates velocity profiles to reflect dual patterns, is the main reason for the outcomes discussed earlier. These significant changes in graphs of principle functions for fractional parameters suggest that fractional models present a more thorough explanation of natural phenomena because the details from the earlier step are captured and used in the system in the following step due to the memory properties of the implemented fractional operator. These results indicate that fractional models offer efficient control over boundary layers however, classical models don't have such features.
To thoroughly analyze the significance of applied boundary conditions, three-dimensional demonstrations of concentration, thermal, and flow fields are provided in Fig. 4a-c respectively. It is clear from Fig. 4b that the thermal function has comparatively lower values if slip effects are taken into consideration. Figure 4c contains two regions corresponding to ramped and constant flow cases at the boundary. In the blue region, the flow profile keeps changing the starting point as long as the value of time ( τ ) changes. This process continues in the domain 0 < τ ≤ 1 . Afterward, the starting point of the flow profile is constant, corresponding to the constant part of the applied flow condition. The respective figure communicates that time variations greatly influence the flow profile for the ramped condition therefore, the use of this condition is helpful for adequately controlling the flow. Figure 5 is constructed to conduct a comparative inspection of flow and thermal profiles for column, lamina, brick, blade, and tetrahedron shapes of nanoparticles. Figure 5a describes that the inclusion of lamina shaped magnesia and graphene nanoparticles in engine oil provides the highest thermal curve. On the other end, the lowest thermal profile is observed when brick shaped particles are considered. The thermal profile www.nature.com/scientificreports/ exhibits the same pattern for slip and no-slip cases. The shape impacts are included in the mathematical model through the shape factor "s", which depends on the sphericity of nanoparticles. The proportion of the sphere's surface area to that of real nanoparticles with equal volumes is known as sphericity. The lamina and blade shaped nanoparticles substantially improve the thermal properties of engine oil therefore, the temperature function of the resulted hybrid nanofluid specifies the highest profiles for these shapes. Oppositely, hybrid nanofluid has a relatively weaker thermal conduction ability for brick and tetrahedron shaped nanoparticles. The velocity fields  www.nature.com/scientificreports/ for constant and ramped cases are respectively displayed in Fig. 5b,c. It is perceived that the order of flow profiles for five distinct shapes is identical to that of thermal profiles. In other words, the flow has the maximum velocity for the suspension of lamina shape particles. This profile is respectively followed by blade, column, tetrahedron, and brick shaped particles. These results express that the distribution of tetrahedron and brick shaped particles signifies the viscous effects. On the opposite end, hybrid nanofluid is less viscous when particles have lamina or blade-like shapes. Hence, it flows with greater velocity as depicted in Fig. 5b,c. www.nature.com/scientificreports/ Figure 6 is created to compare temperature and flow distributions for different combinations of engine oil, magnesia, and graphene nanoparticles. In this figure, Eo-MgO is a magnesia based nanofluid, and relevant graphs are obtained by substituting Υ Gra = 0 in mathematical relations. Likewise, graphene based nanofluid is denoted with Eo-Gra. The figures are prepared for this case by placing Υ MgO = 0 in the final solutions. Eo-Gra-MgO is the main hybrid nanofluid of this work, which contains uniform fractions of both nanoparticles. According to Fig. 6a, the temperature distribution of engine oil receives maximum augmentation when both nanoparticles are dispersed in equal proportions. Contrary to this, pure engine oil shows the lowest values for the thermal function, www.nature.com/scientificreports/ which is somewhat obvious considering its insignificant thermal attributes. In comparison to the thermal conductivity of magnesia, graphene possesses a substantially greater thermal conductivity therefore, the temperature graph of Eo-Gra is relatively higher than that of Eo-MgO. Figure 6b,c reveal that the flow velocity of engine oil is higher than the velocities of other investigated combinations. This highest flow curve is followed by Eo-Gra, Eo-Gra-MgO, and Eo-MgO in the respective order. The conspicuous difference in densities of magnesia, engine oil, and graphene is the primary cause of such velocity patterns. In addition, the immersion of nanoparticles in conventional fluids yields more viscous fluids. Therefore, regular engine oil, due to its lowest density and comparatively weaker viscous nature, indicates the highest flow speed. However, amalgamating engine oil with uniform concentrations of magnesia and graphene considerably disturbs its density. As a result, hybrid nanofluid has a lower density when equated with the density of Eo-MgO. On the contrary, its density is higher than that of graphene based nanofluid. Figure 7 illustrates how the boundary layers of flow and temperature are affected by the rise in total proportion ( Υ ) of magnesia and graphene particles. Figure 7a depicts that enhancements in the total proportion elevate the temperature graph. This figure further describes that the solution of the heat equation has minimum magnitudes for a zero value of Υ , which indicates that nanoparticles have no physical involvement. These notable differences in thermal field graphs demonstrate that the heat-transportation propensity of pure engine oil is highly ineffective for industrial processes. However, when engine oil is hybridized with magnesia and graphene nanoparticles, its heat-carrying capacity is boosted because of the strong intrinsic features of suspended nanoparticles, which improves the functionality of the emerging hybrid nanofluid. The ensuing hybrid nanofluid, because of the boost in thermal features, absorbs heat comparatively faster and in a greater amount; therefore, temperature fluctuations at the boundary wall occur expeditiously. Consequently, Fig. 7a depicts higher thermal profiles that signify a considerable temperature rise. As far as the flow distribution is concerned, the inverse conduct of the velocity profile is discerned from Fig. 7b,c. The corresponding figures further communicate that the immersion of magnesia and graphene nanoparticles leads to producing a pronounced drop in the flow velocity. The viscosity of the host fluid keeps increasing as long as the concentration of nanoparticles continues to rise, which is one of the pivotal characteristics of nanoparticles in the physical sense. This viscosity augmentation results in reduced flow speed therefore, a decline in velocity profile is seen in the respective figures. In addition, temperature slip effects also decelerate the flow. The response of the flow distribution under the relative dominant and weak actions of several forces, like viscous, and thermal and diffusive buoyancy forces, is reported in Fig. 8. Due to the varied characteristics of these forces, flow development is either aided or resisted. In this work, Gr 1 symbolizes the thermal Grashof number, which is influenced by the temperature gradient. In the mathematical sense, the thermal buoyancy force indicates a direct association with Gr 1 , whereas viscous forces and Gr 1 share an inverse correspondence. Similarly, the mass Grashof number is characterized by Gr 2 , which shows a direct connection with the diffusive buoyancy force and is associated with viscous forces in an opposite manner. Figure 8a,b exhibit that positive alterations in Gr 1 heighten the flow graphs. An identical result for flow patterns against the modification of Gr 2 can be perceived from Fig. 8c,d for constant and ramped velocity cases, respectively. Physically, the elevating values of Gr 1 and Gr 2 specify that the boundary wall has a comparatively enhanced temperature and concentration gradient is strong. The conventional currents eventually emerge as a consequence of concentration changes and additional heating disturbing the density. Ultimately, the viscous force is left with a negligible contribution since convectional currents not only yield the buoyant force but also assist in augmenting its intensity. Thus, the deformation faces no significant resistance; therefore, the hybrid nanofluid indicates a greater velocity, and a rise in the corresponding curve can be followed from Fig. 8. A comparative study of concentration, velocity, and heat equations for regular and fractional models is conducted by the dint of Fig. 9. Figure 9a,b describe that thermal and concentration solutions derived using the fractional approach have minimum outputs than those established employing the conventional model. Furthermore, there is no influence of the slip condition on the temperature field in this regard. However, the graphical outcomes of the velocity distribution are very interesting because the respective condition also has a significant contribution in this case. For the uniform condition, the behavior of the velocity field is identical to the aforementioned behaviors of thermal and concentration functions. In this case, the fractional-order solution demonstrates a profile lower than that representing the solution procured via the classical model. Whereas, the ordering of velocity profiles changes when a ramped condition is considered. In this case, the graph of the velocity solution having fractional order is higher. Furthermore, it is discerned that applying the temperature slip condition provides lower graphs of energy and flow distributions as compared to the graphs prepared in the absence of this condition. The results in Fig. 9 also endorse the fact that fractional models, because of their memory features and order-adjustment characteristics, are more effective for an extensive and accurate description of physical mechanisms.
A comparative evaluation of Nusselt number (Nu) for slip and no-slip effects and also for multiple shapes of working hybrid nanoparticles is carried out with the aid of Fig. 10a. It is spotted that improvement in Nu is not the same for each type of considered shape. For instance, the immersion of lamina shaped magnesia and graphene nanoparticles yields the maximum increment in the value of Nu. Contrary to that, the minimum boost in heat transfer rate is witnessed when particles have brick shapes. On the basis of this observation, it is concluded that the most effective shape of nanoparticles to ameliorate the thermal efficacy of industrial fluids is the lamina shape. Furthermore, the consideration of slip effects lowers the outcomes of Nu. Figure 10b,c are provided to investigate Nu and skin friction coefficient ( C f ) for multiple single and dual combinations of magnesia and graphene nanoparticles with engine oil. A comparative inspection reveals that Nu and C f behave differently when the particle's allocation is maximized. Precisely speaking, raising the inputs of Υ Gra and Υ MgO from 0.00 to 0.02 results in the escalation of Nu and diminution of C f . In Fig. 10b,c, the specific proportion of nanoparticles for each presented combination is mentioned along the y−axis. Figure 10b displays that when magnesia and graphene particles have maximal and identical fractions ( Υ Gra = 0.02 = Υ MgO ), Nu produces the highest outcome in comparison to other examined amalgamations. Figure 10c depicts  www.nature.com/scientificreports/ C f is for nanofluid having magnesia nanoparticles. Table 4 is organized to chew over the implications of changing the volume concentrations of hybrid particles for the thermal usability of engine oil. Moreover, the amelioration in heat transport rate is anticipated in terms of percentage. It is found that a slight enhancement in Υ induces a pronounced rise in Nu. Table 4 indicates a boost of 33.37% in thermal effectiveness when magnesia and graphene nanoparticles attain maximum fractions ( Υ Gra = 0.02 and Υ MgO = 0.02 ) during hybridization. This improvement in Nu is quite substantial and supports the employment of hybrid nanofluid that is being analyzed in processes where one of the essential focuses is the efficient cooling of conduits. www.nature.com/scientificreports/ To fully understand the involvement of each considered shape of nanoparticles in strengthening the thermal potential of carrier fluid, Table 5 is developed for seven different values of Υ . These results indicate that laminashaped particles are the most significant when it comes to uplifting the thermal features because Nu is the maximum ( Nu = 8.1363 ) for this shape. In this regard, the performance of blade-shaped particles is comparatively weaker than that of lamina-shaped particles. It is noticed that Nu for brick-like shapes of particles has the lowest outputs as equated to Nu outputs corresponding to other shapes. On making a percentage-based comparison, it is found that the heat transfer rate is only improved up to 8.50% when particles are of brick shapes. These improvements for blade, column, and tetrahedron shapes are 17.79%, 13.91%, and 9.24% in a respective sequence. These noteworthy percentage differences highlight the momentous influences of particles' shapes on maximizing the industrial functionality of ordinary fluids. These findings emphasize that the shapes of embedded hybrid particles are essential factors for the optimization of insufficient thermal characteristics of traditional fluids. Based on the provided results, it can be concluded that the information extracted through theoretic analyses without accounting for shape impacts may not be fully reliable for practical uses. Table 6 is prepared to examine the influence of parameters ζ , σ , and η on Nu for both slip and no-slip cases. It is identified that Nu increases due to the enhancement of these parameters. The outcome is the same whether slip effects are taken into consideration or not. In Table 7, the computational results for Sherwood number (Sh) are compared for two dissimilar inputs of Schmidt number (Sc), taking into consideration the impact of fractional parameters. The table reports that Sh gets escalated when the magnitude of Sc is high. In this case, the diffusion coefficient is small and the viscous impacts are dominant therefore, a boost in the mass transfer rate occurs. Moreover, the parameters ζ , σ , and η tend to enhance the value of Sh. The changes in C f because of these parameters are investigated by the dint of Table 8. This inspection is conducted for ramped and constant boundary velocity functions. The table describes that C f follows inverse patterns for two considered cases against the escalation of ζ . For the constant case, C f gets attenuated, whereas C f produces rising values for the ramped case. This pattern is also followed, subject to increasing alterations of σ . As far as the contribution of the parameter η is concerned, C f communicates declining values for the ramped case as well as the constant case. These enhancements and diminutions of Nu, Sh, and C f in Tables 6, 7, and 8 are purely dependent on the kernel of the applied fractional operator. The results communicated www.nature.com/scientificreports/ in these tables accentuate the fact that the fractional model established in this investigation offers more effective control over heat transfer and flow processes as contrasted to regular models. The specificity and correctness of outcomes can be assured via adapting such models by making the necessary modifications to fractional parameters. The considered problem, which involves heat transfer and flow over a vertical surface, has multiple real-life applications, and the presented findings are useful in this regard. For instance, flow over the vertical surfaces of buildings generates wind loads that can affect the structural design and safety of the building. Regarding heat exchangers, flow over the vertical surface is important for heat transfer between fluids. The flow can enhance the heat transfer rate by promoting mixing between fluids and increasing the surface area available for heat transfer. The flow over the vertical surfaces of wind turbine blades has a substantial contribution to generating lift and producing power. Similarly, the flow over vertical cliffs and shorelines can cause coastal erosion by carrying away sediment and rock. Moreover, flow over a vertical surface is commonly used for cooling electronic devices like computer chips, CPUs, and other electronic circuits. These are a few of the numerous applications of the considered problem. The presented outcomes improve the understanding of flow and heat transfer over a vertical surface. With a better understanding of these phenomena, architects and engineers can design better buildings that are more energy-efficient and comfortable for occupants. This will also reduce energy consumption. These results also help the engineers in the development of such cooling systems that remove heat from the surface more effectively. With a better understanding of heat transfer and flow over a vertical surface that this study offers, the fuel efficiency of ships and automobiles can be enhanced. Moreover, these results suggest that engine oil hybridized with magnesia and graphene particles is a more useful fluid for the lubrication of machines as compared to regular fluids. This is also one of the practical applications of our results. Since our results consist of exact solutions, they can be used to verify the numerical techniques formulated for solving the fractional-order models. Our results also enable the possibility of getting a suitable agreement between theoretical outcomes and the experimental data by using the order-variation property of the involved fractional operator. www.nature.com/scientificreports/

Conclusion
The primary objective of this work is to forecast the escalation in engine oil's thermal capacity because of its hybridization with magnesia and graphene nanoparticles. To thoroughly inspect the involvement of nanoparticles in disturbing the flow patterns and ameliorating the thermal and material characteristics, nanoparticles are assumed to have column, brick, tetrahedron, blade, and lamina shapes. To conduct this analysis, a fractional model is established by the dint of generalized Fourier's and Fick's laws. Ramped flow and temperature slip conditions are jointly taken into account for the first time. The basic governing system is composed of flow, concentration, and energy equations. The Prabhakar fractional operator is implemented to include a multi-parametric  www.nature.com/scientificreports/ kernel into thermal and diffusion flux equations, changing the classical system into a fractional one. The inclusion of dimension-free quantities into basic equations and the application of the Laplace transform are two principal steps for the procurement of solutions. The variations in boundary layers and fall and elevation in profiles of velocity and thermal functions are explicated via graphs. A comparative report on the performance of nanoparticles for several shapes is provided. Moreover, the augmentation in Nusselt number in terms of percentage, shape features, and consequences of modifying the fractional parameters are also investigated. Some comparative illustrations of primary functions, extracted from standard and fractional models, are produced to emphasize the critical impact of fractional techniques for modeling purposes. A brief inspection of shear stress is conducted for fractional parameters and several combinations of graphene, magnesia, and engine oil. The key observations of this analysis are summarized as • Engine oil's hybridization with equal proportions of magnesia and graphene nanoparticles provides a 33% amelioration in its thermal efficiency.  Table 5. Improvement in Nu due to various shapes of graphene and magnesia particles.  www.nature.com/scientificreports/ • When hybrid particles are evenly immersed, their material properties and shape effects escalate the viscosity, due to which the boiling point of hybrid nanofluid rises. Consequently, its potential for heat transportation enhances, and it possesses higher thermal stability. • An increment in collective volume fraction leads to raising the profile of the thermal field. However, the flow profile indicates an inverse trend. • A significant variation in skin friction coefficient for small modifications of fractional parameters demonstrates that the fractional model can adequately control the shear stress. • The lamina shaped hybrid nanoparticles provide the highest values of Nusselt number.
• When fractional parameters are varied, the flow profile demonstrates inverse patterns for constant and ramped cases. • Graphene nanoparticles have a markedly greater influence on strengthening the thermal features than magnesia nanoparticles do. • The buoyancy forces substantially accelerate the flow of hybrid nanofluid.
• In contrast to the slip condition, the thermal curve is higher for the no-slip temperature condition.
• Due to memory features, generalized Fick's and Fourier's laws describe diffusion and thermal fluxes more effectively. • Nusselt number reveals that the highest heat transfer rate is specified by the presented hybrid nanofluid as equated to that of other observed nanofluids and pure engine oil. • The paired employment of the fractional model and ramped velocity function offers improved flow control.

Future research recommendations
• The basic model of this study can be modified to investigate flow problems for other geometries such as disks, cylinders, channels, and pipes. • This model can be extended for two and three-dimensional problems.  www.nature.com/scientificreports/ • Some new results can be obtained by operating other fractional derivatives for the same problem, and comparative analyses can be conducted. • The combinations of other rate-type fluids with different nanoparticles can be studied using the appropriately modified version of this model.